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The discrepancy between maximally entangled states and maximally non-classical quantum cor¬ 
relations is well-known but still not well understood. We aim to investigate the relation between 
quantum correlations and entanglement in a family Bell inequalities with A-settings and d outcomes. 
Using analytical as well as numerical techniques, we derive both maximal quantum violations and 
violations obtained from maximally entangled states. Furthermore, we study the most non-classical 
quantum states in terms of their entanglement entropy for large values of d and many measurement 
settings. Interestingly, we find that the entanglement entropy behaves very differently depending on 
whether N = 2 ox N > 2: when A = 2 the entanglement entropy is a monotone function of d and 
the most non-classical state is far from maximally entangled, whereas when A > 2 the entanglement 
entropy is a non-monotone function of d and converges to that of the maximally entangled state in 
the limit of large d. 
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I. INTRODUCTION 

Quantum correlations arising in space-like separated 
measurement events can be subject to strong correlations 
that cannot be explained by classical physics [I]. The 
strength of classical correlations is bounded by Bell in¬ 
equalities, which are known to have quantum violations, 
certifying the non-classicality of the physics at play. A 
necessary condition to violate a Bell inequality is shar¬ 
ing an entangled quantum state [2]. For pure states, the 
standard measure of entanglement is the entanglement 
entropy [3], defined as the von Neumann entropy of a 
subsystem of the state. The relation between entangle¬ 
ment and the strength of quantum correlations has long 
been considered an interesting question. 

For the simplest Bell scenario, with two parties Al¬ 
ice and Bob performing one of two two-outcome mea¬ 
surements, the CHSH inequality [4] constitutes a tight 
bound on the set of classical correlations [5]. To maxi¬ 
mally violate the CHSH inequality, one must distribute a 
state with maximal entanglement entropy (the so-called 
maximally entangled state). However, when the CHSH 
inequality was generalized to the CGLMP inequality [6] 
which considers scenarios with Alice and Bob choosing 
one of two d-outcome measurements, it was discovered 
that the maximal quantum violation of the CGLMP in¬ 
equality for some low values of d > 2 cannot be achieved 
by distributing a maximally entangled state between Al¬ 
ice and Bob [7, 8]. The CGLMP inequality has been 
studied for scenarios with number of outcomes ranging 
up to the order of 10®, and the most non-classical state 
is found to be far from maximally entangled [9]. Recent 
results [10] suggest that the discrepancy is due to the 
chosen measures quantifying quantum correlations. 

Quantum correlations are relevant in Bell scenarios be¬ 


yond the scope of the GGLMP inequality. Significant re¬ 
search effort has been directed at Bell scenarios where 
Alice and Bob can choose between more than two mea¬ 
surements. Such scenarios with two-outcome measure¬ 
ments for Alice and Bob have been explored in e.g. Refs. 
[11, 12] and many-outcome scenarios in e.g. Refs. [12, 13]. 
In addition. Bell scenarios with many settings and more 
than two parties have been studied in [14]. 

Here we present a much simplified version of the multi¬ 
setting generalization of the CGLMP inequality derived 
in Ref. [15] i.e. for scenarios where Alice and Bob per¬ 
form one of N d-outcome measurements. We investigate 
quantum violations of this inequality for arbitrary N and 
d, using both maximally entangled states and general en¬ 
tangled states. In addition to the maximal violations of 
the inequality, we study the entanglement entropy of the 
most non-classical quantum state by both analytical and 
numerical means. 

Interestingly, we find that the nature of the most non- 
classical state is very different depending on whether N = 
2 (the CGLMP inequality) or A > 2. When A = 2, the 
entanglement entropy of the most non-classical state is 
a monotonically decreasing function of d and thus as we 
increase d the state becomes less entangled. However, 
we find that when A > 2 the entanglement entropy is a 
non-monotone function of d and that the entanglement 
entropy of the most non-classical state converges to that 
of the maximally entangled state in the limit of large d. 

II. MULTI-SETTING AND MANY OUTCOME 
BELL INEQUALITIES 

Let Alice and Bob share a bipartite quantum state of 
local Hilbert space dimension D: |'0) € 0 C'°. Al¬ 

ice and Bob can each perform one of A measurements, 


2 



FIG. 1: A graphical illustration of the chaining associated 
with the inequality (1) where the arrows represent the relation 
“<’ between two measurement outcomes. 


indexed by x^y = 0,...,N — 1 respectively, with associ¬ 
ated outcomes , By each of which can take d different 
values, where d < D, on some ordered space. We will 
focus on cases where the number of outcomes equals the 
local Hilbert space dimension; that is when D — d. The 
constraint D = d will from now on be assumed unless 
otherwise is stated. Scenarios with D > d will be briefly 
discussed. 

For the given Bell scenario, we can construct the fol¬ 
lowing Bell inequality, 


N-l 

BN,d = P{An-1 > Bq) -\- P{An < Bn) + 

n—0 

N-2 

+ Y, PiBn+l < An) > 1 . ( 1 ) 

n—0 

The inequality (1) is a chained Bell inequality, and we 
illustrate it graphically in Figure 1. 

To prove this assertion, we extend the method of Ref. 
[9]. First note the following fact: {Am-i > Bm-i} n 

> A[4_2} n ... n {Bi > Aq} n {Aq > Bq} c 

> Bq}. Therefore, if we take the complement 

of both sides, it holds that {Am-i < Bq} C {A^v-i < 
Bn-i}^{Bn-i < ^7V-2}U. . .U{Bi < Ao}U{Ao < Bq} 
(see Figure 1). Any classical probability distribution 
in the given Bell scenario takes the form of a convex 
mixture of fully deterministic probability distributions, 
for which all probabilities are trivial i.e. either zero or 
one. Therefore, since randomness is not a fundamen¬ 
tal property but a mere statistical property in classical 
probability distributions, it must hold that P{An-i < 
Bq) = 1 — P{An-1 > Bq) < P{Aq < Bq) + P{Bi < 
Aq) P{B]s[-i < An- 2 ) -t- p\Ai<i-i < Htv-i), from 

which the inequality (1) follows. 

Our inequality is (up to a permutation of parties and 
measurement settings) a simplihed version of the inequal¬ 
ity derived in [15], analogous to the simplihcation of the 


CGLMP inequality in Ref. [9]. An advantage of our in¬ 
equality is that its form is independent of d. Another 
advantage is that the inequality only depends on the rel¬ 
ative ordering of the outcomes, and not their ascribed 
values. 

For sake of completeness, we show that of the specihc 
case when the outcomes take the values 0,..., d — 1, our 
inequality is equivalent to that of Ref. [15] (again, up 
to a permutation of parties and measurement settings), 
namely 

lN,d = ([.Bq —Ao])-|-([Ao —i3i])-|-...-l-([RAr-i~^Jv-i]) 
+ {[AN-i-BQ-l])>d-l (2) 

where {X) = J2'i=o IPi^^ = 0 [^] = ^ mod d. 

Simply note that we can write [X] = X — dlX/d\ 
which implies that the inequality can be written ( — 

1 - d{ En=0 liBn - An)/d\ + L(^n - i?„+l)/dj + 

[(Ajv -1 — Bq — l)/dj)) > d — 1. If we use that Ai,Bi e 
{0,..., d — 1}, the inequality (1) is implied. 


III. QUANTUM VIOLATIONS 


We now investigate the quantum properties of the 
Bell expression We consider both violations us¬ 

ing maximally entangled states i.e. a state of the form 
\'4’max) = l/v^Sfc=o and non-maximally entangled 
states. 

Let us denote the outcomes of Alice’s and Bob’s mea¬ 
surements hy a, b G {0,..., d — 1}. In Ref. [15], particular 
measurements were introduced with the motivation that 
the maximally entangled state has the property that if 
Alice measures an observable with an eigenstate |s), and 
Bob measures an observable with conjugated eigenvec¬ 
tors |s)*, correlations will be perfect. Therefore, we let 
Alice and Bob measure in the respective bases given by 




A;=0 

d-1 




/c=0 


( 3 ) 

( 4 ) 


where ai^'^ = x/N and = (1 — 2y)/2N for x,y = 
0,..., N — 1, and oj = . 

We have conducted extensive numerical studies using 
semidefinite programs (SDPs) [16] to investigate whether 
these measurements are optimal. However, the problem 
of optimizing the quantity Rw.d over the measurement 
of Alice, the measurements of Bob, and the shared state, 
cannot be directly maped into an SDP. Instead we have 
used a see-saw technique in which we first use an SDP 
to optimize only over the measurements of Alice. Then 
we run another SDP to optimize over the measurements 
of Bob while fixing Alice’s measurements as obtained in 
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the previous SDP. Finally, we fix Alice’s and Bob’s mea¬ 
surements as returned by the two previous SDPs and run 
a third SDP optimizing over only over the shared state. 
This procedure is iterated until a stable value of is 
obtained. Naturally, this procedure does not guarantee 
that B]\[^d achieves its global optimum under quantum 
theory. Yet, for all investigated cases ranging over sev¬ 
eral low values of N,d, the obtained value of B^^d by 
the see-saw method returns the above measurements of 
Alice and Bob. Due to this numerical evidence, we will 
conjecture the above measurements as optimal. 

The numerics supporting the conjecture of the optimal 
measurements is performed with the constraint that D = 
d. We have used SDPs in an analogous see-saw technique 
to investigate the value of Bff^d in various scenarios when 
D > d. Specifically, we have considered N = 2,3,4, 5, 
d = 2,..., 9 and D = d, ...,10. We have observed no 
instance in which D > d yields an enhancement over the 
case D — d. 

Following the conjecture of the optimal measurements, 
we can define the Bell operator, ^N,d, associated with 
the Bell expression Bw.d as 


^N,d = ^ |a)A.Ar-i(aU,Af-i <8) \b)B,o{b\B,o 

a'>b 

N-1 

+ EE I y4.,n (o-l ^ |^)s,n(^|s,n 

n—0 b>a 
N-2 

+ EE |a)yl,n(a|^,n ® \b) B,n+l{b\B,n+l- (5) 

n—0 a>b 

Thus, we can write that = Tr {p£§N,d) where p is 
some shared quantum state. The optimal value of Sjv.d 
is found from computing the smallest eigenvalue of !^N,di 
and the corresponding eigenvector is the associated quan¬ 
tum state that gives rise to the maximal violation i.e. it 
is the most non-classical quantum state. 

However, in the case of pure quantum states we can 
extensively simplify the eigenvalue computation by using 
the Schmidt decomposition to write any pure quantum 
state in the form |(^) = Y^’l^Q^k\kk) for Afc being real 
numbers respecting the normalization X]fc=o 
is then possible to write the Bell expression Hjv.d in the 
following form 


d-l 

BN,d = ^ MklXkM 

k,l=0 


where the d x d matrix M is given by 


Mki = N6ki 




( 6 ) 


( 7 ) 


This expression is derived in appendix A. 

The problem of finding the maximal quantum violation 
of the inequalities reduces to finding the smallest eigen¬ 
value of M. However, we will first investigate optimal 


values of Sjv.d achievable by performing measurements 
on a shared maximally entangled state. 


A. Violations from maximally entangled states 

The maximally entangled state, \ipmax) has Schmidt 
coefficients = l/Vd for k — 0, ...,d— 1. Inserting 
this into (6), the value of the Bell expression is easily 
computed from 


d-l 


-^A/'.ddV’max)) — , ^ ^ 


k,l—0 


( 8 ) 


Let us study the limiting case of many measurements, i.e. 
find the value of B^oAl^max)) = limjv^oo BN,di\'4’max))- 
We first compute the probability distribution arising 
from Alice and Bob making measurements (3) and (4) 
on \iprnax)- This can be obtained directly from equation 

(Al): 


PAr,d(a, 6|x,y) 


2 sin^ (n (^a-b + ai^'^ + ^ 

sin^ (^a - b + 


(9) 


Using the probability distribution, we can expand the 
Bell expression BN,d{\i^rnax)) in (1) as 




Sin 


2 / TT 


— ) 
2NJ 


d2sin2(-(i_^)) 




2N - 1 


sinM5(i-2v)) sin^ (5 (j + 1- 2^)) 


\2NJ 
■ (10) 


In the limit of large N, it is evident that there is a contri¬ 
bution to Bao^dil'ipmax)) Only from the final term in the 
sum associated to j = d — 1. Using Taylor expansion of 
the sine function, we find that the limit is 

Boo,d{\‘^‘max)) — (H) 

which goes to zero as d —?► oo. In fact, the right-hand-side 
of (11) is the no-signaling bound of B^^d- This is realized 
from the fact that the no-signaling bound of lN,d in (2) 
coincides with its algebraic bound, that is zero, and that 
the affine transformation relating Bj^^d and lN,d is lN,d = 
d X Bf^^d — 1 up to permutations of parties and settings. 
Thus, in the limit of many measurement settings, the 
above result of quantum theory does converge to the no¬ 
signaling bound. However, the algebraic bound of B]\[^d 
is zero, and thus there is a discrepancy between the no¬ 
signaling bound and the algebraic bound of B^^d- 

In the other limiting case of many outcomes, finding 
the value of BN,oo(.\tp7nax)) = limd_>oo ^Ar.dd^mao:)) IS a 
bit more involved. We can use the fact that the sum in 
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equation (8) can be evaluated as an integral, which can 
be related to the Trigamma function, defined as 

the second logarithmic derivative of the Gamma function, 
^ logr(z), as follows 


rl rl sin 

BN,oo{\i’max}) = N-N / - 

Jo Jo S 


^ (N-l)Tr(x-y) ^ 

sin {tt{x — y)) 




2N 


Sin 


/ TT N 

K2NJ 




dxdy 


( 12 ) 


In appendix B, we prove this relation. It is worth noting 
that for the special case of iV = 2, it was shown in [9] 
that the limit evaluates to i? 2 ,oo(|V'maa;)) = 2— 16Cat/7r^ 
where Cat is Catalan’s constant, corresponding to what 
was obtained for the CGLMP inequality in Ref. [6]. 

Evidently, unless we let both N,d —>■ oo, we cannot 
reach the algebraic bound of by using maximally 

entangled states. This puts our inequality in contrast to 
that of Ref. [15] in which the algebraic bound coincides 
with the no-signaling bound and thus can be asymptoti¬ 
cally achieved in the limit N ^ oo. 


B. Optimal quantum violations 

We will now investigate the optimal quantum viola¬ 
tions of the inequalities (1) by sharing general entangled 
states. As we have already mentioned, the optimal quan¬ 
tum violation of the inequality corresponds to the small¬ 
est eigenvalue of the matrix M and the entangled state 
giving rise to the violation is the associated eigenvector 
of M. This reasoning relies on the measurements (3) and 
(4) indeed being optimal choices. To strengthen this con¬ 
jecture, we have verified the optimality of these results 
numerically using SDPs [16] for the cases N,d = 2 ,..., 10. 

In Figure 2, we have computed the maximal violation 
of the inequalities, the violation using a maximally entan¬ 
gled state, and also the violation due to the approximate 
state which is discussed below. We have chosen 

some particular values of = 2, 3,10,100 and d ranging 
from d = 2 up to approximately about d = 0.5 x 10® in 
exponentially increasing increments. The computation 
of the smallest eigenvalue of M (the maximal violation) 
can be efficiently implemented using power iteration over 
the Krylov space associated to M. We elaborate on the 
numerical technique in appendix C. 

The results suggest that the optimal violation of the in¬ 
equality approaches the algebraic bound of namely 

zero, as d —>■ oo for all investigated values of N. The 
convergence to the algebraic bound is more rapid with 
increasing N. 


C. Approximating the most non-classical state 

Due to the apparent discrepancy between maximally 
entangled state and the most non-classical state, J'i/'opt), 



FIG. 2: The optimal quantum violations for N = 2,3, 10,100 
and dimensions up to about d — 0.5 x 10®. Solid lines are 
obtained from maximally entangled states, the dotted lines 
from the most non-classical quantum state, and the dashed 
lines from the approximation state (13)- 


in terms of their ability to minimize B^.d, it is interesting 
to investigate the properties of J'i/’opt)- 

For simplicity, if we represent {ipopt) in the Schmidt 
basis as JV’opt) = J2k=o^k\kk), it follows from the sym¬ 
metries of the matrix M that Itpopt) is subject to the 
symmetry property that Afc = Xd-i-k whenever d > 2. 
However, despite being easy to compute for given val¬ 
ues {N,d), the exact form of Itpopt) is difficult to find by 
analytical means. 

However, we propose an approximation of the state 
\ipopt) for arbitrary N and d, which denote by [‘I’i^p): 



1 \kk) 

v^^KfcTi)(d-fcF^ 


(13) 


where ^ Kd + 1)(^~ j)]~” is the normaliza¬ 

tion. 

The special case amounting to = 2 was initially 
proposed in Ref. [17] to approximate the state maxi- 
mally violating the CGLMP inequality. The state |$app) 
was extensively studied for dimensions up to the order 
of 10® in Ref. [18]. On the one hand, it was shown that 

/o\ 

14>app) yields a good approximation of the optimal viola¬ 
tion of the CGLMP inequality, especially for d < 1000, 
and that it yields a violation which converges to the alge¬ 
braic bound in the limit of many outcomes [18]. On the 
other hand, the numerics of Ref. [18] shows that is 

quite a bad approximation of the state |'0opt) in terms of 
entanglement entropy. Although it was shown that the 
entropy of |4>ipp) becomes F(|$ipp)) = 1/2 in the limit 
of many outcomes, where it converges to the most non- 
classical state [18]. For general N, the numerics in Figure 
2 indicate that the violation due to the approximate state 
|4>i^p) converges to the algebraic bound of Bf^^d in the 
limit of many outcomes. Furthermore, it appears that 
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FIG. 3: The entanglement entropy of the most non-classical 
state (solid line) and the approximate state (dashed 

line) for N = 2,3, 10,100 and dimensions up to d = 0.5 x 10®. 


with higher N, the violation due to increasingly 

well approximates the optimal quantum violation of the 
inequality. 

The question of how well the state |$aS) approximates 
the entanglement entropy of the most non-classical state 
IV'opt) is investigated in the next section. 


IV. MAXIMAL NON-CLASSICALITY VERSUS 
ENTANGLEMENT 

In order to investigate the relation between maximal 
non-classicality and entanglement, we will study the en¬ 
tanglement entropy of the state \ipopt) by extensive nu¬ 
merics for high N and d. We will also study to what 
extent the degree of entanglement of \ijjopt) can be ap¬ 
proximated with the introduced state |d)i^p). 

Using the Schmidt basis representation, the entangle¬ 
ment entropy of a state \(p) is given by 


d-l 

E{\<l))) = (14) 

k=0 

which is easily realized to be a number between zero and 
one (see note [22] for conventions of units). The entangle¬ 
ment entropy is maximized for the maximally entangled 
state for which E{\4>max)) = 1- 

In Figure 3 we plot the entanglement entropy of the 
state \ijjopt) and the approximate state for a few 

values of N and for dimensions up to about d = 0.5 x 10®. 

Remarkably, our results show that for the investigated 
values of N exceeding N = 2, the entanglement en¬ 
tropy of the most non-classical state is not monotonic 
with d. The sharpness of the trough of the entangle¬ 
ment entropy increases with the value of N. Our numer¬ 
ics suggests that E{\tpopt)) is a non-monotonic function 
whenever N > 3. 


(o') 

In addition, while for small values of d the state I'happ) 
provides a less good approximation of the entanglement 
entropy of the most non-classical state [18], we see that 
for N >3, I'bi^p) becomes a significantly better approxi¬ 
mation of the entanglement entropy of IV'opi)- In fact, the 
accuracy of the approximate state with respect to {ipopt) 
appears to increase with both N and d. 

Due to the very slow convergence of the entanglement 
entropy of IV'opt), especially for fix N and variable d, it 
is hard to tell how it behaves for d beyond the range we 
have studied. However, it is reasonable to assume that 
as d —>■ oo, the entanglement entropy converges to that of 
the state |$i^p). In appendix D, we show the following 
limit 

lim E(\^^^'>)) = I if IV = 2 ^ , 

Evidently, the entanglement entropy behaves very dif¬ 
ferently depending on whether we consider case of the 
CGLMP inequality {N = 2) or allow for more measure¬ 
ments. Also, this result implies that the entanglement 
entropy (at the very least for the state |4>i^p)) is a non¬ 
monotonic function of d whenever N > 2, which falls in 
line with our numerical results in Figure 3. 

The above results show that (i) entanglement entropy 
of the optimal and approximate state is non-monotonic 
for V > 2 and that (ii) it converged to that of the max¬ 
imally entangled state for N > 2. Note however, that 
this does not imply that the optimal state or approxi¬ 
mate state are the maximal entangled state for A > 2 
and d —>■ oo. It simply means that the difference between 
their corresponding entanglement entropies becomes in¬ 
finitesimal. An intuitive picture can be obtained by plot¬ 
ting the corresponding Schmidt coefficients Xi as a func¬ 
tion off, as was done in [9]. One observes that the shape 
of the curve is [/-shaped and that for increasing d, as 
well as N the middle part becomes flatter and flatter. 
The above result suggests that as d —>■ oo the function 
X{x) = Xixd] remains [/-shaped over x S [0,1] for N = 2, 
while for A > 2 it becomes essentially a constant func¬ 
tion for X € (0,1), as in the case of the maximally entan¬ 
gled state, apart from a peak with infinitesimal support 
at a; = 0 and x = 1. In fact, those two peaks with in¬ 
finitesimal support make the maximal violations of the 
Bell inequality for this state differ from that of the max¬ 
imally entangled state. 

Note that entanglement entropy as defined in (D6) is 
the commonly used measure for entanglement of a state. 
However, we can also look at the relative entropy or 
Kullback-Leibler divergence between the maximally en¬ 
tangled state and the optimal or approximate state, i.e. 
KL('ipmaxl(/>) where |^) is either the optimal or the ap¬ 
proximate state. This is shown in Figure 4 as a function 
of d for different values of A. Furthermore, it is shown 
in appendix D that 


lim KL(iP^axlHpJ) = 

d—foo 


jlogCjv 


N 


A = 2, 
A > 3, 


(16) 
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FIG. 4: The Kullback-Leibler divergence, for fixed values of 
N, between the maximally entangled state \ipmax} and \(j)}, 
which for solid lines represents the optimal state \ipopt} and 
for dashed lines represents the approximate state 

with Cn = 2w-i0fr(l - 2/N)/T{3/2 - 2/N). Note 
that the divergence for = 2 is very slow, O (log log d). 
This expression agrees with the numerical results shown 
in Figure 4. Furthermore, it is in accordance with the 
above picture, where we can attribute the divergence of 
the relative entropy for iV = 2 to a finite difference of the 
\i and 1/d for infinitely many i, corresponding to a finite 
interval of a; G [0,1], while the finite value of the relative 
entropy for TV > 2 reflects the fact that and 1/d differ 
only for finitely many 1, corresponding to the two peaks 
at a; = 0 and a: = 1 which have infinitesimal support. 

To complete the discussion, we plot 
limd^oo KL('!/max|‘hiS) ^s a function of N in Fig¬ 
ure 5. One sees that the relative entropy between the 
maximally entangled state and the approximate states 
for infinite d goes to zero as TV —>■ oo, which is in 
accordance with the results of [15], meaning that the 
peaks of A (a;) at a; = 0 and x = 1 become smaller as 
N increases until for TV —)■ oo the approximate state 
becomes identical to the maximally entangled state. 


V. CONCLUSIONS 

In this work, we have introduced a much simplified ver¬ 
sion of the mutli-setting generalization of the CGLMP 
inequality first presented in Ref. [15]. We have conjec¬ 
tured the optimal measurements and provided numer¬ 
ics supporting this claim for some low values of (TV, d). 
Then, we studied quantum violations of our inequalities 
using both maximally entangled and generally entangled 
states and demonstrated that the maximal violations are 
obtained by non-maximally entangled states. We gave 
strong numerical evidence that the Tsirelson bound of 
the Bell expression (1) coincides with the algebraic 
bound of that quantity, namely zero. However, our main 
result is the strong evidence for the non-monotonicity of 



Settings N 


FIG. 5: The Kullback-Leibler divergence between the max¬ 
imally entangled state [ipmax) and the approximate state 
14?IS) S'® S' function of the number of the number of settings 
TV in the limit d ^ oo. 


the entanglement entropy of the most non-classical state, 
as a function of d for TV > 2. 

We see that the most non-classical state for any TV > 2 
become more and more entangled with increasing d and 
its entanglement entropy converges to that of the max¬ 
imally entangled state. This is very much opposite to 
the quantum behavior observed for the CGLMP inequal¬ 
ity [18] i.e. when TV = 2. Our results suggest that the 
discrepancy between entanglement and the strength of 
quantum correlations is strongly dependent on the num¬ 
ber of settings of the Bell scenario and is eradicated in 
the limit of large d, given more than two settings. Fi¬ 
nally, recent work [19] studying the ability of random 
pure states to violate the CGLMP inequality has shown 
that with increasing d, the inequality is increasingly vio¬ 
lated on average. Extending such an analysis to the case 
of TV > 2 settings would be of interest. 
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Appendix A: Computing the Bell Operator 

Here we derive Bell operator, M, from equation (7), 
associated to the measurements (3) and (4). 

Given any shared quantum state written in the 
Schmidt basis as ](/) = X)fe=o ^k\kk), the associated prob- 
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ability distribution P]\[^dici,b\x,y) can be simplified to 

k.l—0 

We can use the probability distribution to evaluate the 
Bell expression in equation (1). If we additionally observe 

that +/3i+i = “ 2 ^ and + 

= 1 ~ 2 ^ > one sees that all probabilities appearing 
in (1) are independent of the measurement choices. By 
direct insertion, we write the Bell expression as 


Bna = ^ X! 


k,l—0 


'^{d- 


3=0 


d-1 

(N - 1) J2id - j>('=-'Kj-27v) 
i=i 

j=i 


(A2) 


where the Bell operator M appears in the square brackets 
(scaled by d^). 

In case otk = I, corresponding to the diagonal elements 
of M, it is easily found that 


- 1 

Mkk = N --—. (A3) 

d 

To evaluate the off-diagonal elements, corresponding to 
k—l ^ 0, we first use the geometric series to evaluate each 
of the three series appearing inside the bracket of (A2). 
Combining the exponentials in trigonometric functions 
leads to the final expression 


d-1 

BN,d = ^ AfeA; 
/c,Z—0 


N5ki 


TV sin 

^ (A^—l)7r(fc — 


■ 

^ dN J 

^ sin 

(^) 

. 


(A4) 


where one identifies the matrix M with the term inside 
the bracket. 


Appendix B: Violations from Maximally entangled 
states as d —>■ 00 

Here we compute the limit i?Ar,oo(|'0maa:)) by evaluat¬ 
ing the integral in equation (12). To compute the inte¬ 
gral, we use the Dirichlet series for the cosecant function, 
namely 


1 

sin (z) 


i=i 


(Bl) 


If we define z = 7r(a; — y), we obtain an expansion 
for one of the two factors appearing in the integrand 


of equation (12) in terms of the variable z. Addition¬ 
ally inferring the notation fi = {N — l)7r(x — y)/N and 
/2 = (2j — l)7r(x — y), we can write the integral expres¬ 
sion in (12) as a sum over an infinite series of integrals 
over integrands of the form sin (/i) . Using Euler’s 

formula, we can write each integral in the infinite se¬ 
ries as a sum of four elementary integrals which can eas¬ 
ily be solved; each over a term of the form cos (/i i f 2 ) 
or sin (/i ± / 2 ). However, two of the four integrals will 
vanish: it is a straightforward computation to find that 
fo fo ^ •^ 2 ) dxdy = 0 for all values of j. 

The task of computing is reduced 

to evaluating the following series of integrals: 


N-N 1 fo Jo (cos (/i - / 2 ) - cos (/i -b / 2 )) dxdy 
The two elementary integrals appearing in this expres¬ 
sion are evaluated to 

^ 27V3(1-cos(|))) , , 

N /cos (/i ± / 2 ) dxdy = \ (B2) 

Jo Jo TT^ 1 t 27 +A) 


where we have defined = j and j_ = j — 1. Using this 
expression, one has 


Bn, 00 — N + 'y ( 
i=i 


2jV3(l_cos(^)) 
7r2 (1 - 2jNf 


2V3(i_cos(^)) 
^^7r2(l + 2(j-l)iV)^' 


(B3) 


The first summation, which we denote Si, can be re¬ 
written as Si = Vsin^ (^) /tt^ J^JLi (j “ 2 lv)”^- 
call that the Trigamma function, ipi(z), can be defined 
through the series ipiiz) = + j)~^) which is pre¬ 

cisely the form of the summation in Si. If we additionally 
use the recurrence relation ipii^ + 1) = i^ii^) ~ Bi 
can be simplified as 



(B4) 


A similar analysis can be done for the second sum in (B3) , 
now denoted S 2 - First using the recurrence relation, and 
then using the reflection property of the Trigamma func¬ 
tion, ipi{z)+ipi{l — z) = -K^ j sin^ (ttz), allows one to write 
S^ = N-Nl3.'^^iii^{YN)il,i{\-^). 

Plugging the final form of Si and S 2 back into (B3), 
yields 


B 


N,oo — 




(B5) 


which is the final expression (12). 


Appendix C: The numerical technique for finding 
maximal violations 

In this paper we have numerically computed both the 
maximal violation of the inequality (1) and the entangle¬ 
ment entropy of the associated most non-classical state 

















for values of d ranging up to about d = 0.5 x 10®, for 
a few values of N. This was achieved by using similar 
techniques as in [9] which we will elaborate on here. 

As previously mentioned, the maximal violation of 
the inequality (1) is found from computing the small¬ 
est eigenvalue of the matrix M (7), and the most non- 
classical state is the associated eigenvector A of M. Re¬ 
call that the components of this eigenvector A are just 
the Schmidt coefficients. Since M is a d x d matrix, the 
computational complexity of the problem scales as 0{d^) 
with increasing d, becoming less and less feasible. Nev¬ 
ertheless, there exist more efficient numerical techniques 
for solving the problem, especially since we are only inter¬ 
ested in the smallest eigenvalue of M and the associated 
eigenvector. 

We have used the method of power iteration over a 
Krylov subspace which is the essence of the Arnoldi and 
Lanczos algorithm [20, 21]. The Krylov subspace of order 
r -|- 1 associated to some n x n matrix A and some n x 1 
vector b is the space spanned by the image of b under 
the first r powers of A. For large r, these images of b 
approach the eigenvector of A associated to its largest 
eigenvalue (in modulus). 

For our purposes, we have fixed r = 20 and started 
the iteration from the Schmidt coefficients associated 
with the maximally entangled state i.e. A*^®^ = l/y/d x 
(1,..., 1)^. However, since we are interested in the small¬ 
est eigenvalue of M, we have decomposed M as 

M = Nld^d - M' (Cl) 

where /dxd is the d-dimensional identity matrix. We now 
exploit the simple fact that the eigenvectors of M are the 
same as the eigenvectors of M' and that the correspond¬ 
ing eigenvalues are related by m = A7 — m', in particular 
for m being the smallest eigenvalue of M and m' being 
the largest eigenvalue of M'. 

The above described algorithm is summarized as fol¬ 
lows: 

1; procedure Power iteration for optimal state 

2 : initialize: A^®i ^ l/v^ x (1,...,!)^ 

3: for k in 1,.., 20 do 

4: 

5: end for 

6: normalize: ^ 

7: H = Af-(A(20))'^M'A(20) 

9: end procedure 

Importantly, to avoid memory problems, we only store 
the current vector A^^^ in memory when running the it¬ 
eration. In addition, since M has a banded structure we 
can store its independent elements in a d/2-dimensional 
vector. Further optimization can be achieved by exploit¬ 
ing the symmetry of M or M'. The program is imple¬ 
mented in C-I--I- and ran on the Arcus cluster at the Ad¬ 
vanced Research Computing, University of Oxford. Data 
and code are available on request. 


Appendix D: Entanglement entropy of the 
approximate state I'f'iS) 

In this appendix we derive the limit (15). That is, we 
evaluate 


1 


d-l 


if” I I (Dl) 


where the Schmidt coefficients are given by 




(D2) 


and = J2j=o [0 + 1)('^ ~ j)] ^ is the normalization. 

We first compute the normalisation to leading order by 
approximating the sum as an integral. For convenience, 
we introduce e = l/d and the rescaled variable x = j/d. 
Then the normalization can be written as 




dx 


'o [(cc-I-e)(l - x)]^ 


(D3) 


where the dots refer to next-to-leading-order terms. Ex¬ 
plicit evaluation of the integral yields to leading order 




2e| logej -f ... 

+ ... 


AT = 2, 
AT > 3, 


where 


4 1 

Cn = 2N-1 




r(i- 2 /A) 

r(3/2-2/A) 


(D4) 


(D5) 


which is finite for A > 3. Note that the explicit form of 
Cn is not relevant for this derivation, but will be relevant 
later. 

Approximating the sum in the entanglement entropy 
in a similar manner one arrives at 


^(l„w)) = M££!f) 

I logej 

2 A/N-i 


pl-e 


+ 


_ log[(x + e)(l-x)] 

A^jlogejio [(a; + e)(l_a;)]w 


dx .... (D6) 


To take the limit d —>■ oo, or equivalently e —>■ 0, of 
(D6), we calculate the leading order contribution from 
both terms. Using (D4) one obtains for the first term 


log (^e-4/^) 

I log e| 


(D7) 


The integral in the second term can be evaluated as 
follows for A = 2 


l-e 


log [(x-f e)(l - x)] 
[(x-f e)(l-x)]^ 


dx = -| logep -I- ... 


(D8) 
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while for TV > 3 one has 


where the number of outcomes goes to infinity. That is 


log.i) 

Jo [(x + e)(l — x)] ™ 


(D9) 


Combining this with (D4), we see that the second term 
in (D6) becomes 


lim = i™ 

a—>-oo d—^oo 



d-1 \ 

^logdAl 

t=o / 


(D12) 


Following a similar strategy as above one obtains 


2 ^ ' log [(x + e)(l - x)] 

A^^llogelio [(,c + e)(l-x)]^ 

+ - N = 2, 
\ TV > 3, 


(DIO) 


/ 4 \ 2 

= + — J log[(x + e)(l-x)]dx + ... 

= log(AAe^-w) + A(-2 + 2e-2eloge) + ... (D13) 


which for TV > 3 becomes zero in the limit e —>■ 0. Plug¬ 
ging (D7) and (DIO) into (D6) and taking the limit finally 
gives 


hmiTdc&W)) 


TV = 2, 
TV > 3, 


(Dll) 


which completes the derivation. 

For completeness we also calculate the relative entropy 
or Kullback-Leibler divergence between the maximally 
entangled state and the approximate state in the limit 


In the limit where d —)■ oo, or equivalently e —> 0, the 
second term converges to —4/TV. Furthermore, we see 
from (D4) that the first term diverges for TV = 2, while 
it converges to logC^v for TV > 3. In total one has 



oo 

logCjv 


N 


TV = 2, 
TV > 3, 


(DI4) 


where is given in (D5). Note that the divergence for 
TV = 2 is very slow. Indeed from the above expression, 
one sees that it diverges as O(loglogd) when oo. 
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